Topological features of proteins from amino acid residue networks 
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Topological properties of native folds are obtained from statistical analysis of 160 low homology 
proteins covering the four structural classes. This is done analysing one, two and three- vertex joint 
distribution of quantities related to the corresponding network of amino acid residues. Emphasis on 
the amino acid residue hydrophobicity leads to the definition of their center of mass as vertices in 
this contact network model with interactions represented by edges. The network analysis helps us 
to interpret experimental results such as hydrophobic scales and fraction of buried accessible surface 
area in terms of the network connectivity. To explore the vertex type dependent correlations, we 
build a network of hydrophobic and polar vertices. This procedure presents the wiring diagram 
of the topological structure of globular proteins leading to the following attachment probabilities 
between hydrophobic- hydrophobic 0.424(5), hydrophobic-polar 0.419(2) and polar-polar 0.157(3) 
residues. 

PACS numbers: 89.75.-k, 87.15.Aa 
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I. INTRODUCTION 

It has long been recognized that limitations on mod- 
eling proteins are due to the complexity of the intra- 
molecular interactions and the lack of an adequate de- 
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scription of the interactions with the solvent |1J. The 
complexity of the interactions has led to sophisticated 
numerical simulations of all-atom models. Thus, it is 
primordial to consider the basic features of the pro- 
tein dynamics to obtain a reliable and less computation- 
ally demanding working model to understand the folding 
process and protein stability. Many theoretical models 
have been proposed to the protein structure prediction 
with reduced degrees of freedom Q. This leads to a 
coarse-grained approach where all-atom interactions are 
replaced by centers of force. 

Beyond the dynamical approach, many natural and ar- 
tificial interacting systems have been modeled as a graph, 
where vertices represent the system elements and edges 
represent their interactions. In particular, some graphs 
present scale free behavior and are called complex net- 
works 0, IH Q Theoretical studies have focused on the 
understanding of their large-scale organization [|| 0] and 
growing mechanism possibly driven by evolutionary 
principles 0, EH ED • 

In this view, networks of interacting amino acid pairs 
have been constructed given the Cartesian coordinates 
of the main chain a-carbon atoms as vertices. Edges are 
established when any other C a atom is within a cutoff 



Euclidean distance R r 
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taken to define a contact network of amino acids. 
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Another approach to model an amino acid network can 
be devised if one considers the center of mass of their side 
chains as vertices. Those centroids are obtained with 
the heavy atoms spatial coordinates of the amino acid 
residues and act as interaction sites. Again, a contact 
network is defined if a centroid is within an Euclidean 
distance R c = 8.5 A of another one and this case is stud- 
ied in this paper. This approach emphasizes the amino 
acid residue hydrophobicity, which is the main force lead- 
ing to a folded conformation 0, ^| . 

With this network modeling, following the work by 
Greene and Higman |l3j. we are able to gather quan- 
titative results characterizing the native conformation 
topology of globular proteins. As we show, this proce- 
dure helps us to rationalize the formation of native states 
from a topological point of view. This is achieved with 
a network description for the inter- residue interactions, 
which helps to interpret experimental results such as hy- 
drophobic scales and the fraction of buried accessible sur- 
face area. Moreover, we evaluate the attachment proba- 
bilities between hydrophobic-hydrophobic, hydrophobic- 
polar and polar-polar residues from this simple network 
model. 

The definitions of the topological quantities used in 
this paper are presented in Sec. [n] with the numerical 
results followed by a discussion. The final conclusions 
are presented in Sec. IIIII Appendix lAl complements our 
study through the distribution of inter-residue contacts. 
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II. RESULTS AND DISCUSSION 



To explore the network model of globular proteins, sta- 
tistical analysis have been carried out for networks built 
from 160 high-resolution and low homolo gy s tructures, 
presenting different folds and chain lengths [2(| , collected 
in the Protein Data Bank (PDB). Chains with residues 
omitted in the PDB database have not been included 
in our data set. This avoids bias in the calculation of 
topological quantities. Our data consist of four sets of 40 
different native-state conformations representing the four 
broad structural classes, all-a, all-/3, a + (3 and a//3 ac- 
cording to the SCOP classification |2jJ. The selection of 
protein structures has been made to permit a numerical 
analysis of the dependence of network topological prop- 
erties on secondary structure contents and their network 
size N, i.e., the total number of amino acids. 

Figure ^ shows an illustrative example of the network 
structure with the vertices calculated from side chain cen- 
troids for the protein 1E0L, collected in our all-/3 data 
set. 



A. One vertex analysis: mean network degree and 
hydrophobicity scale 

Figure [3 shows the average vertex (degree) connectiv- 
ity (k) for the four structural classes as a function of N. 
It is obtained as an average over the connectivities fcj, 
where fej is the number of edges connecting the vertex i to 
other vertices in the network. The value of (k) increases 
quickly for small proteins and shows an asymptotic trend 
for all classes. This figure indicates a predominant lower 
average connectivity for all-a conformations compared 
with all other classes. 

The typical (k) for proteins are to be compared with 
other biological systems Q : S. cerevisiae protein-protein 
network, (k) = 2.39 and N = 1870; E. coli metabolic 
network, (k) = 7.4 and N — 778; neural network of C. 
elegans, (k) — 14 and N — 282; Silwood Park food web, 
(k) = 4.75 and N = 154; Ythan estuary food web, (k) = 
8.7 and N = 134. 

A protein mapping onto a contact network clearly in- 
troduces constraints on the network topological proper- 
ties. The first one is due to the finite size of real data, 
which introduces a cutoff degree depending on N [22| . 
Second, size and steric considerations introduce a more 
relevant constraint in the topology of amino acid net- 
works, for example, limiting the possible values for the 
degree k. Thus, the degree dependence must be charac- 
terized by a limiting finite value for the maximum degree 
kmax, irrespective of the chain size. This can be observed 
for our data sets in the inset of Figure |21 The hitting of 
the maximum degree of a vertex can be translated as 
performance of how to connect vertices in the folding 
process. Once this value is reached, the network rewires 
edges connected to that vertex or stops this specific ver- 
tex dependent process. 
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FIG. 1: Example of the amino acid residue network for the 
PDB code 1E0L, a small 37-amino acid protein in the all-/3 
class. Figure (a) shows the ribbon diagram obtained with 
RASMOL visualization tool. Figure (b) displays the corre- 
sponding network with their numbered vertices representing 
the center of mass of the side chains. This figure was drawn 
using VISONE (www.visone.de). 



This conclusion is also reinforced by a similar study 
carried out by Vendruscolo et al. 01 with the "betwee- 
ness" Bi, which correlates with kf, when comparing na- 
tive and Monte Carlo generated transition states. As ex- 
plicitly demonstrated for six proteins, highly connected 
residues in the native state may not correspond to the 
highly connected ones (key residues) in the transition 
states. 
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FIG. 2: Average connectivity as a function of the network 
size N for the structural classes: all-a (o), all-/3 (□), a + [3 
(+) and a/(3 (A). Inset: the maximum connectivity k ma x 
for each network as a function of N. 



Protein structure is sequence dependent. Thus, differ- 
ent of other real network models, where vertices are in 
general considered to be of the same type we are dealing 
with systems, which result from the intrinsic properties 
of the amino acid residues. 

To obtain a deeper insight of the role played by each 
residue in forming native conformations, we investigate 
how to relate their character to the ability of making 
contacts with other residues. For this end, we display 
in Fig. |21 the average connectivity for each residue one- 
letter code from all 160 structures, taking into account 
the number of times each residue appears in all sequences. 

The sequence of amino acid residues displayed in Fig- 
ure |3] was chosen because it groups their characteristics: 
aliphatic non-hydrogen-bonding (Gly, Ala, Pro, Val, He, 
Leu), aromatic non- hydrogen-bonding (Phe), sulfur con- 
taining (Met, Cys) and hydrogen-bonding (Trp, Tyr, His, 
Thr, Ser, Gin, Asn, Glu, Asp, Lys, Arg) (23J. 

Figure^ shows a clear trend for lower average connec- 
tivity played by amino acid residues in proteins classified 
in all-a structural class. Nevertheless, the highest av- 
erage connectivity is presented by residues in the a//3 
protein class. This class has mixed a-helices and par- 
allel /3-strands. Intermediary pattern of connectivity is 
found for residues in all-/? and a + (i protein classes with 
predominant higher connectivity for the f3 class. There- 
fore, it is likely that the higher average connectivity for 
the residues in a/ (3 class also comes from interactions 
among residues in different secondary structures. The 
a + f3 class presents regions of cv-helices and anti-parallel 
/3-strands largely segregated along the chain. We also 
conclude that the connectivity in this protein class re- 
sults mainly from contacts among the elements in each 
secondary structure type. Figure [Sb shows the average 
connectivity of each residue calculated from all networks 
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FIG. 3: (a) Average connectivity of each amino acid residue 
according to the structural class, (b) Average connectivity 
and standard deviation of each amino acid residue calculated 
from all networks. The four more connected residues are Val 
(V), He (I), Leu (L) and Phe (F), while the less connected 
ones are Gin (Q), Glu (E), Asp (D), Lys (K) and Arg (R). 



considered in our data set, stressing their general ability 
of making network contacts. 

Can we relate this ability to any intrinsic amino acid 
property? In the following, we argue that those sharp 
evaluated average values can be understood in terms of 
the concept of hydrophobicity. 

We recall that the hydrophobic effect concept as the 
leading force in folding globular proteins comes from ex- 
perimental evidences. This is illustrated by the find- 
ing of mostly non-polar residues in the core of protein 
structures 0, 0] as a consequence of avoiding the ex- 
posure to water of amino acids with hydrophobic side 
chains. Therefore, the hydrophobic/hydrophilic charac- 
ter of amino acids has been viewed as an important is- 
sue in understanding the formation of native structures. 
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For example, the spatial distribution of hydrophobic and 
hydrophilic residues has led to the definition of score 
functions based on hydrophobic profiles as a mean to 
detect native protein folding among many decoy struc- 
tures 

In spite of many scales ranking the hydrophobic char- 
acter, due mainly to the particular measurement method, 
there is a consensus about the more hydrophobic (Phe, 
He, Leu, Val) and more hydrophilic (Lys, Arg, Asp, Glu) 
residues. Therefore, one expects to obtain more (less) 
connected vertices as the ones ranked in the top (bot- 
tom) of the hydrophobic scales. 

To investigate how the average network connectivity 
(shown in Figure^) is related to the hydrophobic char- 
acter of residues, we have checked its correlation with the 
127 scales collected in Ref. [26J. The most significant cor- 
relations are observed with the well known scales in the 
literature: Kyte and Doolittle, Juretic et al. and Janin- 
2, all presenting correlation coefficient p = 0.94; Taylor 
and Thornton (p = 0.93) and Sereda et al.-2 (p = 0.92), 
following the notation for these scales as in Table I of 
Ref. [26( . On the other hand, non-significant correlations 
are obtained with the scales: Hopp (p = 0.10); Wilce 
et oZ.-3 (p = 0.14); Michel et al., Colonna-Cesari and 
Sander (p = 0.21). 

Although we can find a very good agreement between 
the amino acid hydrophobic caracter and its average con- 
nectivity, Figure[3f a) clearly shows that this quantitative 
character is related to the secondary structure contents. 
Yet, this has experimental consequences on the defini- 
tion of the free energy transfer. Therefore, any method 
used to predict elements of secondary structure based on 
typical hydrophobicity profiles need to be reconsidered. 
Those quantitative values are expected to be correlated 
to the exposure pattern of the amino acids to the aqueous 
environment. This can be readly investigated considering 
only the average connectivity and the available data on 
individual amino acid accesible surface area (independent 
of its location in a peptide). 

Our results show that the average connectivity reveals 
a weak correlation with the average buried accessible sur- 
face area (ASA) A" - {A) (p = 0.39) Hill, where A ° is 
ASA calculation based on Gly-X-Gly models related to 
the unfolded state and (^4) is an average ASA calculated 
from native protein structures. However, the correlation 
coefficient increases to 0.90 for the fraction of buried ASA 
1 — (A) /A , which reflects the relative loss of their acces- 
sible surface areas during the folding process. 



B. Two vertices analysis: assortativity 

So far we have focused the discussion from a single ver- 
tex point of viewing. However, m any real-world networks 
also exhibit degree correlations [2!3,|3(j, i.e., the condi- 
tional probability P(k'\k) that an edge leaving from a 
vertex, with degree k, and arriving at a vertex, with de- 
gree fc', is dependent on k. 



This two- vertices correlation is expressed by k nn (k) = 
Y^ k , k' P(k'\k) and furnishes the average connectivity 
of the nearest neighbors of vertices with connectivity 
Concerning this aspect, the networks are 
classified as to show assortative mixing, if the degree cor- 
relation is positive, a preference for high-degree vertices 
to attach to other high-degree vertices, or disassortative 
mixing, otherwise. Therefore, whenever k nn (k) increases 
with k, one obtains a correlated network, thus showing 
an assortative mixing. On the other hand, a decreasing 
behavior as a function of k shows a disassortative mixing. 

Our numerical analysis of the two-vertices correlation 
shows that the dependency of k nn (k) on the degree k can 
be described by a linear relation for all our networks (data 
not shown), k nn {k) = a + bk, with a, b > 0. Moreover, 
this result leads to the average of k nn (k) over all vertices, 
{knn)N — J2k P(k)k n n{k)j where P(k) is the observed 
degree distribution of each network, to have also a linear 
behavior as a function of (k) because (k) = J2k kP(k). 

We depict in Figure ^ the behavior of (fc„„)jv as a 
function of N and of the structural classes. We observe 
again a stationary value for large N, reflecting what 
we can call the main property of amino acid residue 
network. In fact, that homogeneous network behav- 
ior ((fe nra )jv ~ (k)) is supported by the inset of Fig- 
ure 0] The linear regression for the 160 proteins leads 
to (fc n „)jv = 0.65(7) + 1.012(7)(fc), with correlation co- 
efficient p = 0.99. Therefore, there is a strict correspon- 
dence, exhibited by the above equation, between the way 
{knnjN and the average degree (k) grow with N. As a 
matter of fact, this can also be noted a posteriori by 
the similarity between the positions of the data points 
in Figures [21 and QJ This reveals a specific pattern of at- 

ta.rtimpnts amnnff vertices for all network si^es and it is 
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FIG. 4: Average of k„„(k) over all vertices of degree k, 
(fcnn)iv, as a function of the network size N for the struc- 
tural classes: all-Q (o), all-/3 (□), a + (5 (+) and a//3 (A). 
Inset: {k nn )M as a function of the average connectivity (k) 
for all structural classes. 
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We recall that, although we observe a linear increasing 
of k nn (k) with k, this behavior is only valid up to k ~ 
k m ax due to physical constraints of the networks. 



C. Hydrophobic and hydrophylic vertices and their 
correlation 

So far, we have reinforced the observed process to- 
ward some high-connected vertices as the final result 
of the protein folding within its biological time. How- 
ever protein folding corresponds to a specific mecha- 
nism of wiring in contrast to the collapse of homopoly- 
mers, likely due to the key residues in according to the 
nucleation-condensation model as a highly cooperative 
process 

To explore further the role played by the connectivity 
of residues leading to the assortative mixing behavior we 
re-define our networks as being composed of two types 
of residues, denoted by H and P. H residues are chosen 
as the more connected ones in Figure and P residues 
otherwise, in such way to obtain two data sets of 11 and 9 
residues, roughly corresponding to the three more corre- 
lated scales we have identified in Sec. Ill Al This makes a 
set for H type residues including Thr as being the lowest 
hydrophobic and another one for P type residues start- 
ing with Ser. Differently from bipartite networks, where 
edges connect only vertices of unlike type, our new net- 
works have edges connecting vertices irrespective of their 
type. Now, to proceed the investigation on how the as- 
sortativity arises in native structures we apply the ideas 
of community forming by Newman [3^ |. In this sense, 
we search a better comprehension of the assortativity 
as consequence of a mechanism for network formation, 
which results from preferential attachments among spe- 
cific residues. Firstly, we quantify the observed assorta- 
tive mixing calculating the fraction of edges eij between 
vertices of type i and j — {H,P}). Table I presents 
the averaged values with their standard deviations for 
the fraction of edges in each structural class and the fi- 
nal averages from all data, with chp — epH- It is clear 
the preference for attachments between H — H and H — P 
vertices. Secondly, we calculate the assortativity coeffi- 
cient [2^, r, a correlation coefficient which lies in the 
range [0, 1]. This coefficient hits the value 1 for perfect 
assortative mixing, i.e. when every edge connects vertices 
of the same type, and the minimum value for a random 
mixing. The low numerical values for r in Table I tak- 
ing into account their standard deviations confirms that 
there is almost no assortative mixing by type, i.e. H — H 
and P — P attachments. Hence, the network surrounding 
residue of type H is likely to be another H or P, with 
probabilities presented in Table I. We may interpret those 
sharp probabilities as the result of a compositional equi- 
librium of H and P residues in globular proteins. With 
our hydrophobic scale defined in terms of (k) we calculate 
how dominant the presence of H residues is. For all-a 
class we obtain 0.549 ± 0.007; all-/?, 0.569 ± 0.008; a + 0, 



0.559 ± 0.007 and a//3, 0.579 ± 0.006. 

Although the tendency of H residues to connect among 
themselves and with P residues is a well known fact, 
this approach presents a manner to quantify the pref- 
erential attachment associated to globular proteins. Fur- 
thermore, the tuning of energy couplings in minimal 
model simulations may profit from this residue attach- 
ment quantification in order to be more realistic in de- 
scribing the folding process. 

From a graph point of view, attachment probabilities 
reveal the wiring diagram when including different kinds 
of vertices as a mechanism for network formation. Those 
probabilities also discard the simple picture of a predom- 
inant hydrophobic core. Our results help to form the 
image that an H vertex brings along other H vertices as 
well as other P vertices in forming the native structure. 

Table I also shows the fraction of edges for all 
classes as a function of three typical ranges of protein 
sizes N. It is interesting to note that when one restricts 
the analysis for chain sizes in those ranges, the fractions 
eij in each range are still comparable within the standard 
deviations. This confirms the robust pattern of network 
formation irrespective of their size. 

D. Three vertices analysis: non-random pattern 

To support the above specific attachment among 
{H, P} residues, we have evaluated the clustering coeffi- 
cient (C) over all vertices. While the average connectivity 
gives information on how connected a vertex is, the clus- 
tering coefficient shows the network transitivity, i.e., the 
probability that two nearest neighbors of a given vertex 
are nearest neighbors themselves forming a triangle 
with this vertex. Notice that, in this aspect, the cluster- 
ing coefficient is a measure of three- vertices correlation. 
The maximum value of (C) is 1 in a fully connected net- 
work, it assumes the value {k)/N for a random network 
and it is constant for regular lattices (3j| . Our results for 
(C) as a function of N for the structural classes exclude 
any random character in its formation (See Figure 01 as 
already pointed out 0, without including any size 
dependence analysis. 

Next, we look at the distribution of inter-residue con- 
tacts from a network perspective. However, since the 
main results of this distribution are already known in 
the literature, we have left them to the appendix. We 
include this appendix to show that network modelling 
distinguishes molecular chains with secondary structures 
from polymers and to present a comprehensive study of 
the topological properties of globular proteins. 

III. CONCLUSION 

Different from the usual network analysis dedicated to 
study the properties of a single system, we try to infer the 
main topological features of a class of systems (globular 
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all-a 0.405(9) 0.422(5) 0.173(7) 0.099(8) 

all-/3 0.420(10) 0.426(5) 0.154(6) 0.073(8) 

a + /3 0.428(9) 0.417(5) 0.155(5) 0.091(6) 

a//3 0.443(7) 0.412(4) 0.145(4) 0.091(6) 

all classes and all sizes 0.424(5) 0.419(2) 0.157(3) 0.089(4) 

all classes with TV < 100 0.419(18) 0.424(9) 0.157(11) 0.070(13) 

all classes with 100 < TV < 200 0.438(9) 0.410(5) 0.152(7) 0.096(8) 

all classes with TV > 200 0.420(5) 0.422(2) 0.159(3) 0.091(4) 



TABLE I: Fraction of edges between a vertex of type i (H or P) and a vertex of type j (H or P). Results are presented 
for globular proteins according to the structural classes and for all classes according to the number of amino acids TV. The 
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FIG. 5: Clustering coefficient calculated with our network 
model as a function of TV. 



proteins) characterized by an unknown mechanism which 
drive them to specific tridimensional configurations, the 
native structures. 

Topological parameters have been obtained to charac- 
terize native structures and their organization pattern as 
a result of the underlying dynamical process of forma- 
tion. Our main results permit the following theoretical 
conclusions: (i) the average number of edges leaving a 
given center of force presents a direct relation with the 
hydrophobic character of the residues and to their frac- 
tion of buried accessible surface area; (ii) the way this 
average number of edges (center of forces) grows with 
the protein chain size is also reproduced in the formation 
of inter-residue interactions among the nearest neighbors 
of that center of force; (iii) the interactions among spe- 
cific residues exhibit a well defined pattern described by 
the results presented in Table I. The above patterns re- 
main, irrespective of the protein size, and this may shed 
light on the protein evolution. 
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APPENDIX A: DISTRIBUTION OF 
INTER-RESIDUE CONTACTS 

There is still another important aspect to be consid- 
ered: the attachments distribution in our networks com- 
pared to a full connected one, as driven by bonded and 
non-bonded forces from inter-residue interactions. This 
can be translated as a topological pattern of inter-residue 
contacts categorized into short-, medium- and long-range 
interactions. Short-range interaction contributions corre- 
spond to residues within a distance of two residues from 
each other, medium range to the ones within a distance 
of 3 or 4 residues and long-range contribution is defined 
to come from more than 4 residues away 36] . Figure 
presents the average relative number of observed inter- 
residue contacts as a function of its sequence separation 
\i — j\ between residues i and j along the main chain. It 
is defined as the ratio of the observed number of contacts 
at distance \i — j\ and its maximum number for the full 
connected network, i.e. with TV(TV — l)/2 edges. This 
figure shows that the main contribution comes from the 
first four neighbors, i.e., from short- and medium-range 
contacts, reflecting the chain connectivity property. In 
particular, for all-a chains, we see that from all possible 
contacts at distance 1, only 90.2% occur; at distance 2, 
43.9%; at distance 3, 60.6%; at distance 4, 52.2%; drop- 
ping to 11.8% and 7.7% at distances 5 and 6, respectively. 
Comparing all classes, the short-range contacts present 
similar relative numbers. Large relative number (but less 
than 1) of contacts at distance 1 means that the chain 
connectivity between successive residues (i, i + 1) is not 
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FIG. 6: Average relative number of inter-residue contacts as 
a function of the sequence separation \i — j\. Results are 
obtained for the four structural classes of proteins with N > 
200. 



We observe a preference for medium-range contacts by 
all-a structures and a very low preference by all-/3. The 
highest relative number of contacts for all-a at distance 3 
reveals the pattern of helices characterized by the (i, i+3) 
contacts. Classes a + j3 and a/ (3 present similar aver- 
age relative number of observed inter-residue contacts. 
It is clear the preference of long-range contacts for all-/? 
structures, mainly compared with all-a structures. This 
overall behavior reveals how vertices related to secondary 
structures interact in a network of amino acids as a func- 
tion of the sequence separation. 
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